Breakdown of Luttinger liquid state in one-dimensional frustrated spinless fermion 

model 
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Haldane hypothesis about the universahty of Luttinger 
liquid (LL) behavior in conducting one-dimensional (ID) 
fermion systems is checked numerically for spinless fermion 
model with next-nearest-neighbor interactions. It is shown 
that for large enough interactions the ground state can be 
gapless (metallic) due to frustrations but not be LL. The ex- 
ponents of correlation functions for this unusual conducting 
state are found numerically by finite-size method. 

PACS: 71.10.Pm, 71.10.Fd 



One-dimensional (ID) Fermi systems have a number of 
peculiarities which distinguish them drastically from 3D 
ones (for review, see |l]). In particular, gapless (metal- 
lic) ID systems of interacting fermions never behave as 
normal Fermi-liquid. Haldane [^ have proposed another 
class of universality Luttinger liquid (LL) state, namely. 
It is characterized by the existence of three branches of 
low- lying Bose excitations, density, current, and charge 
excitations with the velocities vs, vj and vn, correspond- 
ingly. The first one is connected with the variation of the 
total energy of the system E under the variation of the 
total momentum P, vs = SE/6P; the second one (f j), 
with the variation of the energy under the shift of all 
the particles in momenta space, that can be done phys- 
ically by the application of magnetic flux to the system 
closed as a ring; and the third one, with the variation 
of the chemical potential /i at the change of the total 
number of particles N, vn — (L/n) Sfi/SN, with L being 
the length of the system. In LL state there are an exact 
relation between the velocities 



X = vjvn/vI = 1, 



(1) 



which is the criterion of LL. The only dimensionless pa- 
rameter which determines all the infrared properties of 
the system (e.g., time and space asymptotics of fermionic 
Green functions and susceptibilities) is the ratio 



vn/vs = vs/vj 



(2) 



Original arguments by Haldane were based on exact 
Bcthe Ansatz solutions as well as on the perturbation 
theory for weakly interacting systems. There are no gen- 
eral and rigorous proof of this assumption but all the 
known analytical and numerical results about ID fermion 
systems confirm Haldane hypothesis pM]. Anderson 0] 
proposed that some 2D systems such as copper-oxide su- 
perconductors also belong to the class of LL which made 



the concept of Luttinger liquid one of the most "fash- 
ionable" in contemporary many-particle physics. There- 
fore the investigation of the status of Haldane hypothesis 
seems to be of importance. Here we present a counterex- 
ample to this hypothesis basing on exact numerical re- 
sults for spinless fermion model. 
We proceed with the Hamiltonian 



H = 



-^E 



c\c,+i 



'-J+l^- 



L L 



=1 



=1 



(3) 



where c\ , Ci are Fermi creation and annihilation operators 
on site i,ni = c\ci. Phase diagram of this model has been 
investigated by us in H. In particular, it has been shown 
that for half-occupied case, p = N/L = 1/2 and arbi- 
trarily small t the ground state turns out to be gapless 
(metallic) along the line V — 2V'. It is the consequence 
of frustrations which lead to the macroscopically large 
degeneracy (finite entropy per site) of the ground state 
at the Ising limit t — 0. A similar result has been ob- 
tained also in Ref. H. It is important that, according to 
our calculations, the metallic region has non-zero width 
in (V, V') plane. One can say with certainty that the 
gap is zero at {V/2) - OM <V' < {V/2). One can show 
also with certainty that the ground state is insulating at 
\{V/2) — V'\ > 1. To check the Haldane hypothesis we 
restrict ourselves by the consideration of the straight line 
V = 2V' where the system is definitely metallic. Simi- 
larly, we have the metallic state for p = 2/3 at V^' = 
and arbitrarily large V or, vice versa, at V^ = and arbi- 
trarily large V . There are rare examples of metallic state 
with strong interactions and it seems to be interesting to 
check the Haldane assumptions for this unusual case. As 
it was already mentioned the original Haldane hypothesis 
is based, on the one hand, on the consideration of exactly 
integrable systems and, on the other hand, on the per- 
turbation treatment of systems with weak correlations. 
Therefore, its validity in the case under consideration is 
not obvious. Note that "unusual" character of metallic 
state at p = 1/2, V « 2V' has been mentioned in Q but 
without specification what this state is. 

We have carried out the calculations of the ground 
state of the model ^ by Lanczos method for finite clus- 
ters with the consequent extrapolation to L — > oo (for 
details, see [p|). Velocities of low- lying excitations has 
been calculated as [Eljll 



vj^^[Ea{L,N)~E^[L,N)] (4) 

ZTT 

VN = - [Eo (L, A^ + 1) - 2£;o (L, iV) + £^0 (i, N - 1)] 

Here Eq {L, N) is the ground state energy of the clus- 
ter with L sites for periodic boundary conditions and N 
particles, Ea{L,N) is ground state energy for antiperi- 
odic boundary conditions (transition to the antiperiodic 
conditions corresponds to magnetic flux $ = 1/2 of the 
flux quantum), Eip{L,N) is the ground state energy for 
minimal nonzero total momentum P = 2Tr/L. Then we 
have verifled the criterion of LL x = 1 using Eq. (yj) . 

The results of the testing calculations for the case 
p = 1/2, V ^ 0,0 < V < 2t where the system has to 
be LL H] are shown in Fig. 1 (open circles and triangles) . 
We also present in the same figure the calculated val- 
ues of X along the line V — 2V' . One can see that at 
V < lot we have, within the accuracy of the computa- 
tions, X « 1, in an agreement with Haldane hypothesis. 
However, for V > 30t the values of x is definitely less 
than unity, that is obvious even without extrapolation 
to L ^ oo since x(^) < 1 ^or finite L and diminishes 
with L increase. Therefore we have demonstrated that 
there are one-dimensional conducting systems of inter- 
acting fermions which are not LL. The breakdown of LL 
picture is caused by the competitions of nearest-neighbor 
and next-nearest-neighbor interactions (i.e. frustrations) 
which allow the system to be metallic in the limit of 
strong interactions. A schematic phase diagram is shown 
in Fig. 2. The question is still open whether the transi- 
tion from insulating state to non-LL conducting state is 
the direct one or there exists intermediate conducting LL 
phase. At the same time, our calculations demonstrate 
that for p = 2/3 the relation (||) takes place with the ac- 
curacy of calculations for any values of parameters under 
consideration even along the lines y = or T^' = 0. It 
would be very interesting to understand analytically the 
reason for the difference between these two cases with 
strong frustrations. 

We also have calculated the static correlation functions 



(4,^0) 



G(i?) 

K (R) = {6nB.Sno) 



(5) 



where brackets mean the averaging over the ground state, 
6ni = rii — p. In LL the following asymptotics have to be 
valid at i? > 1 il 



G{R) ^ Y. C^sm[{2m + l)kFR]R-''" 

oo 

KiR) ^ J2 Dm cos {2mkFR)R-'^'" 

m=0 



(6) 



where r?™ = ^-'^'^ + 2 {m+ ^f €'^'^,9^ = 2m^e'^'fi 
(m > 0) ykp = p/2 is the Fermi momentum. The most 
important exponent a determins the behavior of one- 
particle distribution function n (fc) near the Fermi sur- 
face 



n (k) w n (kp) — Csign {k — kp) \k — kf 



(7) 



where a = ryo — 1. 

However we cannot use these expressions a priori since 
for the model under consideration the system is not al- 
ways LL. We have found the asymptotics of the correla- 
tion functions by direct computation. It is known (see, 
e.g., 0) that it is very difficult to find the correlation 
exponents from the calculations for a given L, even as 
large as L = 32. Therefore we use the finite size scaling 
technique pi| . Specifically we use the following proce- 
dure. 

Our aim is to find the function ip{R) = {(j){0)(t){R)) ^o 
for the infinite chain. Direct calculations give us the func- 
tions /(i?, L) = {(t>{0)(l){R))L for R< L. From the sym- 
metry considerations we have f{R, L) = f{L — R, L). Let 
us introduce the function r{R, L) to have, by definition, 
ip[r{R,L)] = f{R,L). Therefore 



lim r{R,L) 

L— i-oo 



R. 



(8) 



Then we introduce the new variable A = R/L so that 
r{R, L) — L ■ r'{X,L) where r'{\,L) is a new unknown 
function. To provide (^ we have limi^oo '''(A,£) = A. 
Also the function r' satisfies the condition r'{X,L) = 
r'(l — X,L). For small A one has r'{X) « A. To satisfy 
all these requirements we try the function r' as a Fourier 
series 



r'(A) 



sin(7rA) + 03 sin(37rA) 4- 05 sin(57rA) 
7r(H- 03 + as -I- . . .) 



(9) 



Using the asymptotic expression similar to (g) for the de- 
pendence /(i?, L) — ip (Lr'(A))at finite L and optimizing 
the result with respect to both On and the correlation ex- 
ponents we can find the latter with high enough accuracy. 
At least the results for the exponents appeared to be ac- 
curate enough for the clusters with 14 < L < 26 used in 
our calculations. For the testing case V — 0,0 <V <2t 
where the system is definitely LL the results for the cor- 
relation exponents coincide with that from Haldane for- 
mula (pi) with the accuracy of 0.5% for the function G (R) 
and 8% for the function K (R). 

In the most interesting case x 7^ 1 we cannot use the 
expression (|6|) and have to restrict ourselves only by the 
consideration of the leading terms in the asymptotics of 
correlation functions which are tried in the following form 



GiR) -(Ci+C2sin(fii'))/i?T 
K{R) ^{D,+D2{-l)'')/R' 



(10) 



(we consider the case p = 1/2). To diminish the number 
of states in the Gilbert space under consideration we use 
only the states which has the same (minimal) energy for 
V = 2V' and i = which corresponds to the considera- 
tion of the case V/t ^ oo,V' /t ^ oo,V/V' = 2.lt allows 
us to consider as large clusters as L = 32. The results of 
the calculations for the correlation functions arc shown 
in Figs. 3, 4. We have found by the technique described 
above 7 = 2.009 ^ 2.013 and 5 = 1.80 ^ 1.83. Note that 
the envelope of the function K{R) turns out to be non- 
monotonous in non-LL regime (see the black circles for 
i? = 2, 4, 6 in Fig.4. 

The results of computer simulation demonstrating pos- 
sible violation of Haldane hypothesis seem to be rather 
unexpected. In particular, we cannot see any simple 
causes for the difference between two frustrated cases: 
p = 1/2, V = 2V' -^ 00 (non-LL behavior) and p = 
2/3, y = 0, y ^ 00 (LL behavior). It would be very im- 
portant to understand these numerical results by regular 
field-theoretical methods. 



schematically by zigzags. 

Fig. 3. The dependence of the correlation functions 
G(R) (Eq.(5)) for L=32; open circles correspond to 
V = V = 0, black ones correspond to y = 2V' , y ^ 00. 

Fig.4. The same as in Fig.3, for K{R) (Eq.(5)). 
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CAPTIONS TO FIGURES 

Fig.l. The dependence of the ratio x (Ec[-(1)) on the 
inverse size of the cluster; empty symbols correspond to 

V = (circles: V = 0.5t, triangles: V ^ IM); black 
symbols correspond to F = 2V' (circles: V — lOt, 
squares: V = 30i, triangles: V = 50i, diamonds: 

V = loot, hexagons: V = 200t). 

Fig. 2. Phase diagram of the model. The bound- 
ary between conducting LL and non-LL phases is shown 
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